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ABSTRACT 

We use high-resolution hydrodynamic re-simulations to investigate the properties of the 
thermal Sunyaev-Zel'dovich (SZ) effect from galaxy clusters. We compare results obtained 
using different physical models for the intracluster medium (ICM), and show how they mod- 
ify the SZ emission in terms of cluster profiles and scaling relations. We also produce realistic 
mock observations to verify whether the results from hydrodynamic simulations can be con- 
firmed. We find that SZ profiles depend marginally on the modelled physical processes, while 
they exhibit a strong dependence on cluster mass. The central and total SZ emission strongly 
correlate with the cluster X-ray luminosity and temperature. The logarithmic slopes of these 
scaling relations differ from the self-similar predictions by less than 0.2; the normalization 
of the relations is lower for simulations including radiative cooling. The observational test 
suggests that SZ cluster profiles are unlikely to be able to probe the ICM physics. The total 
SZ decrement appears to be an observable much more robust than the central intensity, and 
we suggest using the former to investigate scaling relations. 

Key words: cosmic microwave background - galaxies: clusters: general - hydrodynamics - 
methods: numerical 
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1 INTRODUCTION 

Galaxy clusters are ideal probes for studies of large scale structures. 
They have typical masses of order lO^"* — W^^AIq ; within a radius 
of a few Mpc they contain hundreds of galaxies orbiting in a gravi- 
tational potential well due primarily to dark matter (DM). They are 
also filled with ICM: hot ionized gas, typically at a temperature of 
1-15 keV, that can be observed via its bremsstrahlung emission in 
the soft X-ray band, and via the Sunyaev-Zel'dovich (SZ) effect in 
the millimetric band. 

Several observational evidences demonstrate that the picture 
of an ICM in hydrostatic equilibrium in the DM potential well 
of the cluster is an over-simplification: above all, scaling relations 
between various cluster properties, like mass, temperature and X- 
ray luminosity, differ from the self-similar predic t ions (see, e.g. , 
the recent reviews by iRosati et al]|2002l : F\^l2005l : iBorganill 20061 
and the references therein). In order to explain the ICM properties, 
many models have been developed so far, which consider the effect 
of radiative cooling and of non-gravitational heating (mainly feed- 
back from supemovae and AGNs). As the complexity in the phys- 
ical description of the ICM increased, hydrodynamic simulations 
have also improved to include such processes, and have become in- 



valuable tools to compare theoretical predictions with current data 
and to forecast the performance of future experiments. 

On the observational side, important improvements are ex- 
pected for the detection of the SZ effect in clusters thanks to an up- 
coming newgeneration of suitable instruments, like AmQ SPTq 
ACtE SZ^I AMiB^E APEX0, and PlanckQ. In order to cor- 
rectly extract and utilize the information contained in these data it 
will be mandatory to keep under control both the dependence of the 
SZ signal on the ICM physics, and the possible biases induced on 
the intrinsic signal by the details of the observational process. 

In the present work we investigate how different choices for 
the physical modelling of the ICM affect the SZ radial profiles and 
the scaling relations between the SZ flux and other intrinsic proper- 
ties of clusters. To this end, we use a set of high-resolution hydro- 
dynamical simulations, where clusters were re-simulated including 
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different physical processes, from radiative cooling, to star forma- 
tion, energy feedback from supernovae, and thermal conduction. 

After this analysis, we adopt the instrumental characteristics 
of the AMI interferometer and simulate the observational process 
for a few clusters in our sample: this enable us to discuss possible 
biases introduced in our results by a realistic observation process. 

The present paper is organized as follows. After briefly intro- 
ducing in Sect. |2] the SZ effect and the main related quantities, in 
Sect.[3]we describe the general characteristics of the hydrodynamic 
simulations used in our study. In Sect. |4] we obtain fitting formulas 
for the cluster SZ profiles, while in Sect.[5]we discuss the effects of 
various physical processes on the scaling relations between some 
global properties commonly used to describe galaxy clusters. In 
Sect. |6] we perform simulated observations assuming the charac- 
teristics of the AMI instrument, and discuss the robustness of our 
results. Finally in Sect.|7]we draw our conclusions. 



2 THE SUNYAEV-ZEL'DOVICH EFFECT 

The thermal SZ effect dSunvaev & Zeldovich|[l972[) is a distorsion 
in the cosmic microwave background (CMB) spectrum due to in- 
verse Compton scattering of CM B photons by hot ionized gas parti- 
cles (for recent reviews see, e.g.. lBirkinshawll999l : ICarlstrom et al] 
[2002; Rephaeli et al. 2005). Its intensity is given by the Comp- 
tonization parameter y, defined as: 

y = I n^T^dl , (1) 

where ks is the Boltzmann constant, or is the Thomson cross sec- 
tion, c is the speed of light, and rie, Te and me are the electronic 
number density, temperature and rest mass, respectively. As evident 
from Eq.l[T), the y-parameter is directly proportional to the cluster 
pressure P integrated along the line of sight. The change in the 
CMB intensity at the frequency v corresponding to a certain value 
of the Comptonization parameter is: 

A/cMB = ^J^^9(^) y . (2) 

where h is the Planck constant. The frequency dependence is given 
by the spectral shape g{x): 

where x = hu /kBTcms- 

In galaxy clusters the Comptonization parameter y typically 
reaches values of order 10^'' on arc-minute scales, correspond- 
ing to fluctuations in the CMB temperature of few ^K. As a 
consequence, observing galaxy clusters through the SZ effect re- 
quires extremely good sensitivity, in order to obtain sufficiently 
high signal-to-noise ratios. The required sensitivity and sky cov- 
erage will be achieved in the near future, thanks to new dedicated 
instruments. Nevertheless, the resolution of these instruments will 
be not enough to detect the wealth of details which can be seen in 
numerical hydrodynamical simulations: substructures, shocks, etc. 
For this reason, in order to be meaningful and helpful, all analyses 
of the SZ effect based on simulations must deal with the actual in- 
strumental capabilities: this is one of the main goals of this paper, 
as discussed below. 



Table 1. Different physical processes included in our re-simulations 



simulation name 


involved physical processes 


ovisc 


non-radiative gas-dynamics with standard viscosity 


Ivisc 


like ovisc, but with low viscosity 


csf 


cooling, star formation and feedback 


csfc 


like csf plus thermal conduction 



3 PROPERTIES OF THE SIMULATED CLUSTERS 

The hydrodynamic simulations used in the present work were 
carried out using GADGET-2 iSDringel 2005), a new version of 
the paralle l Tree-l-SPH (Smoot hed Particle Hydrodynamics) code 
GADGET l lSpringel et al.ll200ll) . GADGET-2 includes an entr opy- 
conserving formulation of SPH dSpringel & Hemauistll2003) , ra- 
diative cooling, heating by an UV background, and a treatment of 
star formation and feedback from galactic winds powered by super- 
nova explosions. The latter is based on a subresolution model for 
the multi-phase structure of t he interstellar medium, as described in 
ISpringel & Hernguistl ( l2003h . The SPH implementation of thermal 
conduction in GADGET-2, which is both stable and manifestly con- 
serves thermal energy - even wh en individual and adap tive time- 
steps are used - was described in I Jubelgas et alj ( l2004h . This im- 
plementation assumes an isotropic effective conductivity parame- 
terized as a fixed fraction of the Spitzer rate. It also accounts for 
saturation, which can become relevant for low-density gas. 

3.1 The sample of re-simulated clusters 

The set of simulated galaxy clusters used in the following analysis 
includes 1 1 objects, which we divide in two samples: seven low- 
mass clusters, with virial mass A/vir 10"/i*^A'/q (LM sam- 
ple), and four high-mass systems with A/vir ~ (1-3 — 2.3) x 
lO^'^h-'^MQ (HM sample). The cluster regions were extracted 
from a (DM-only) cosmological simulation (see lYoshida et al.l 
l200lh with box-size 479 Mpc of a flat ACDM model with mat- 
ter density parameter flom ~ 0-3, cosmological constant density 
parameter S7a ~ 0.7, dimensionless Hubble parameter h — 0.7, 
power spectrum normalization ag — 0.9 and baryon density pa- 
rameter il h = 0.04. Using the "Zoomed Initial Conditions" (ZIC) 
technique jTormen et al .111 9971) ■ these regions were re-simulated at 
higher mass and force resolution by populating their Lagrangian 
volume in the initial conditions with a larger number of particles, 
while appropriately adding high-frequency modes. At the same 
time, the large-scale tidal field of the cosmological environment 
is correctly described by using low-resolution particles. The un- 
perturbed particle distribution (before displacements) is of glass- 
type (,Whitei,1996.) . Gas was added only in the high-resolution re- 
gion by splitting each original particle into a gas and a DM one. 
Thereby, the gas and the DM particles were displaced by half 
the original mean interparticle distance, so that the centre-of-mass 
and the momentum were conserved. Each gas particle has mass 
JTigas = 1.7 X 10* Mq, while the mass of the DM particles is 
mum ~ 1.13x 10^ A/q ; thus, the total number of particles in- 
side our simulated clusters is between 2 x 10^ and 4 x 10*, depend- 
ing on their final mass. For all simulations, the gravitational soften- 
ing length was kept fixed at e = 30 kpc comoving (Plummer- 
equivalent) for z > 5, and was switched to a physical softening 
length of e — 'ah~^ kpc at 2 = 5. 

Each object of the cluster set was re-simulated four times, in- 
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Table 2. The sample of simulated clusters. For each object we report name, virial mass Afviri virial radius Hvir, mass-weighted temperature T, X-ray 
luminosity Lx in the [0.1-10 keV] band, and central value j/o "f the Comptonization parameter. Clusters are divided in two samples depending on their mass: 
low-mass (LM) and high-mass (HM). 



LM sample 


A^vir [Mq/K] 


-Rvir [kpc//l] 


T [keV] 


Lx [erg/s] 


I/O 




g676 


1.1 X IQi-* 


983 


1.1 


3.2 X 10** 


3.6 X 10- 


-5 


g3344 


1.1 X 10^-* 


1002 


1.4 


2.2 X 10** 


2.9 X 10- 


-5 


g6212 


1.1 X IQi-* 


1000 


1.3 


3.0 X 10** 


4.4 X 10- 


-5 


gl542 


1.1 X 10^* 


982 


1.3 


3.0 X 10** 


4.2 X 10- 


-5 


gib 


4.5 X IQi-* 


1585 


2.6 


4.2 X 10** 


5.0 X 10- 


-5 


glc 


1.9 X lO^-* 


1181 


1.4 


3.0 X 10** 


4.3 X 10- 


-5 


gld 


1.5 X IQi-* 


1098 


1.2 


1.0 X 10** 


1.4 X 10- 


-5 


HM sample 


Mvi, [Mq/K] 


iJvir [kpc/h] 


T [keV] 


Lx [erg/s] 


yo 




g8a 


2.4 X IQls 


2758 


10.6 


8.1 X 10*^ 


7.4 X 10- 


-4 


gla 


1.5 X IQis 


2358 


7.0 


4.8 X 10*5 


4.3 X 10- 


-4 


gSla 


1.3 X IQis 


2276 


6.2 


3.2 X 10*5 


3.5 X 10- 


-4 


g72a 


1.4 X 10^^ 


2299 


5.8 


1.8 X 10*5 


1.9 X 10- 


-4 



eluding a different set of physical processes for the gas component 
(see the list in Table[TJ. Both simulations referred as ovisc and Ivisc 
follow non-radiative gas-dynamics: for ovisc a standard viscosity 
scheme is adopted, while Ivisc runs were car ried out using a mod- 
ified a rtificial viscosity scheme suggested bv lMorris & MonaghanI 
l ll997l) . where each particle evolves in time its individual viscos- 
ity parameter. Whereas in this scheme shocks are as well captured 
as in the standard scheme, regions with no shocks do not suffer 
for residual non-vanishing artificial viscosity. Therefore, turbulence 
driven by fluid instabilities can be much better resolved by the new 
scheme and, as a result of this, galaxy clusters can build up a suf- 
ficient level of turbulence-powered instabilities along the surfaces 
of the lar ge-scale velocity structure present in cosmic structure for- 
mation (iDolag et al. I l2005l) . In simulations which include star for- 
mation and feedback {csf and csfc), the supemovae efficiency in 
powering galactic winds was set to 50 per cent, which corresponds 
to a wind speed of about 340 km s^^. For the csfc simulations, 
which also include the effect of thermal conduction, we assume a 
fixed fraction 1/3 of the Spitzer rate. A more detailed description of 
the properties of simul ated galaxy cluster s including thermal con- 
duction can be found in lDolag et al.l 12004). 

Table|2]summarizes the main properties of the 1 1 clusters, ob- 
tained using the ovisc simulation. In particular, for each object we 
list its virial mass Mvir, virial radius 7?vir, mass-weighted tempera- 
ture T, X-ray luminosity Lx in the [0.1-10 keV] band, and central 
value yo of the Comptonization parameter. Notice that both T and 
Lx are computed inside i?vir, which was defined using the over- 
density thresho ld dictated by the spherical top-hat model (see, e.g., 
lEke et ai] |l996h . The same properties in simulations with a differ- 
ent ICM modelling show discrepancies at most of few per cent for 
virial masses and mass- weighted temperatures. Conversely, typical 
X-ray luminosities, when compared to the ovisc runs, are reduced 
on average by a factor of 2 in the Ivisc simulations and by a factor 
of 5 in both the csf and cfsc runs. 



4 FITTING RADIAL PROFILES 

In this section we discuss the procedure adopted to find fitting re- 
lations for the radial profiles of the SZ maps of our clusters. To 
produce two-dimensional (2D) SZ maps we followed a standard 
procedure which employs the SPH kernel to distribute the desired 



quantities (in this case the Comptonizatio n parameter y), on a grid. 
Details on the method may be found in iDolag et alj j2006b and 
iRoncarelli et al.l ^2006 ). In particular, we considered the simula- 
tion outputs at z = 0, and for each cluster we obtained three SZ 
maps, corresponding to projections along the cartesian axes. The 
SZ maps, 4 IVIpc on a side, were integrated over 8 Mpc along the 
line of sight: we checked that such integration is enough to enclose 
most of the cluster signal, so that no artificial cutoff is introduced. 

Examples of SZ maps are shown on the right column of Fig[T] 
where we present two systems, gla (upper panels) and glc (bot- 
tom panels), representative of the HIVI and LM samples, respec- 
tively. In both cases we show the output of the non-radiative simula- 
tion ovisc. For completeness we also display the corresponding 2D 
maps for the gas density (left column) and for the mass-weighted 
temperature (central column). 

The radial profiles obtained from the SZ maps are strongly 
dependent on the cluster mass: for this reason we decided to obtain 
fitting profiles for the LM and HM samples separately. In order to 
obtain best-fitting relations for the profiles of the Comptonization 
parameter y, we started from three-dimensional (3D) data, and built 
a model for the SZ emission directly linked to the 3D cluster shape. 
In particular, we obtained the 3D pressure profile for each cluster 
of our sample and fit it using a simple analytical formula, which 
we then integrated along the line of sight. As we are going to show, 
the resulting expression, evaluated at the correct integration limits, 
provides a good fit for the SZ radial profiles. 

The 3D pressure profiles of each cluster were obtained by cal- 
culating the mean value of gT for particles in spherical shells of 
equal logarithmic width out to 47?vir. In order to extract a mean 
profile for each mass bin, and for each set-up of the physical pro- 
cesses, we need to appropriately normalize each profile, obtaining 
dimensionless quantities. Assuming a self-similar model for the 
density profiles, differences in pressure are due to differences in 
temperature only. We found that the resulting pressure profiles for 
our simulated clusters, however, do not properly superimpose when 
normalized by their temperature, indicating that also density pro- 
files show non-negligible deviations from a self-similar behaviour. 
For this reason we preferred to use an empirical normalization, 
adopting the value of the pressure integrated up to a given fixed ra- 
dius. We chose i?5oo, the radius at which the mean overdensity of 
the cluster is 500 times the critical density of the universe (roughly 
corresponding to 60 per cent of i?vir)- The normalization factor ob- 
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density temperature SZ 



gla 




-5 -4 -J -2 2*68 10 -7 -6 -5 -4 

log [g/cm^] [keV] log [Compton y] 




-e -5 -* -J 0.5 1.0 15 2-0 -8 -7 -6 -5 



log [g/cm^] [keV] log [Compton y] 

Figure 1. Projected maps of gas density (left column), mass-weighted temperature (central column) and Comptonization parameter y (right column) for two 
galaxy clusters: gla (top panels) and glc (bottom panels), both taken at z = and considering the non-radiative simulation ovisc. The displayed region is a 
square with side 4 Mpc. The colour scale is shown at the bottom of each panel. 



tained in this way for each cluster was also used to obtain the mean 
SZ profiles. 

After obtaining a mean 3D profile for each cluster subsam- 
ple, we looked for a fitting relation that is also integrable along the 
line of sight. We chose the following expression for the 3D fit in 
spherical coordinates [R, 6, 0): 



F{R) = hD{R + rpy + K 



(4) 



where R = y + + is the 3D clustercentric distance (in 
units of the virial radius Rvir) and I3D is the intensity normaliza- 
tion factor. The free parameters used for the fit are three: rp (also 
given in units of Rvir), b and K. 

In order to obtain the SZ map fitting relation we rewrite Eq.Q 
in cylindrical coordinates (r, 9, z) and integrate along z, obtaining: 



f{r) 



F{r, z)dz ■ 



1 2D 



[b-i) 



Kz 



(5) 



where r = + y'^ is the 2D clustercentric distance (again in 

units of the virial radius 7?vir) and h, I2 are the integration lim- 



its adopted. From a theoretical point of view, integration along the 
line of sight should be performed from the centre to infinity. For 
practical reasons, however, the outer limit can be safely set to a ra- 
dius above which the intensity does not significantly contribute to 
the result: in our case we choose I2 = 2 x Rvh- As for the inner 
limit, the 3D profile in the simulations can be robustly computed 
only down to radii enclosing at least of order of 100 particles, so 
our fitting profile for i? — » is actually an extrapolation. In partic- 
ular, the contribution from the cluster centre to the integral along 
the 2-axis leads to errors in the SZ profile. However, if we give up 
integrating down to 2; ^ 0, we can minimize the contribution from 
R and get a correct integrated profile. For this reason we used 
h = 0.015i?vir and we accounted for the flux underestimation in 
the factor J215. 

As an example, in Fig.[2]we show radial profiles for the pres- 
sure P and for the Comptonization parameter y (left and right 
panels, respectively). Results are displayed for the physical model 
ovisc and for both cluster samples: LM (upper panels) and HM 
(lower panels). Solid lines indicate the mean of the logarithmic pro- 
files, with corresponding la error bars, while dashed lines show the 
fitting relations obtained using Eqs.lO and l|5j. We notice that the 
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Figure 2. The radial profiles for 3D pressure P (left panels) and for the Comptonization parameter y (right panels) as obtained from the ovisc simulations. 
Upper and lower panels refer to the LM and HM samples, respectively. Solid and dashed lines are the mean of the logarithmic profiles and our fitting relations, 
respectively. Error bars show error on the mean. The smaller box below each panel shows the mean of the logarithmic residuals with the corresponding root 
mean square (rms) value. 



3D fit is not always within the error bars: this is mainly due to the 
noisy tail of the pressure profile outside i?vir- Our interpretation 
is confirmed by looking at the logarithmic residuals shown in the 
smaller box below each panel: here we plot their mean and rms val- 
ues. At any rate, we will see that SZ profiles are more regular, as a 
result of the integration process, and the integrated relation will be 
an extremely good fit. 

Extending our analysis to the simulations which include dif- 
ferent physical processes, we computed the fitting parameters 
which let us derive the 3D fitting profiles according to Eq.l|4j, and 
the 2D ones according to Eq.l|5j- The results are summarized in 
Tabled 

The resulting SZ profiles for all physics are plotted in FigH] 



for both mass bins. We find that differences due to cluster physics 
are mainly visible in the low-mass clusters, and affect in particular 
the inner region of the profiles: therefore such differences are un- 
likely to be probed through observations (see the discussion below). 
We also note that, for both mass bins, the models csf and csfc have 
the same rp, b, and K, and only slightly differ in intensity J215. 
This is due to the fact that thermal conduction acts mainly locally, 
so it does not significantly affect quantities averaged on sufficiently 
large scales, like the radial SZ profile. 

As we stated above, there is a significant difference between 
the mean profiles obtained for high- and low-mass clusters, as 
clearly shown in Fig. [3] We note that there is also a selection effect 
which may be responsible for this effect: the HM clusters of our 
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Table 3. Fitting parameters for 3D and 2D profiles (Eqs.(4) and (5), respectively). 



Cluster sample 


Physics 


/3D[erg/cm3] 


-^215 




rp 


b 


K 


LM sample 


ovisc 


7.48 X lO"^"* 


1.83 X 10" 


7 


0.045 


-2.75 


0.0 




Ivisc 


8.93 X 10-1"* 


1.88 X 10" 


7 


0.075 


-2.9 


0.0 




csf 


1.00 X 10-13 


2.36 X 10" 


7 


0.06 


-2.55 


0.0 




csfc 


8.10 X 10"" 


1.99 X 10" 


7 


0.06 


-2.55 


0.0 


HM sample 


ovisc 


6.37 X 10-13 


3.94 X 10" 


6 


0.05 


-2.6 


-0.45 




Ivisc 


6.56 X 10-13 


5.26 X 10" 


6 


0.08 


-2.6 


-0.45 




csf 


6.56 X 10-13 


3.52 X 10" 


6 


0.09 


-2.8 


-0.45 




csfc 


6.43 X 10-13 


3.74 X 10" 


6 


0.09 


-2.8 


-0.45 




10"^ 



HM sample 




.00 



Figure 3. Comparison between the mean logarithmic SZ profiles obtained 
from cluster simulations including different ICM modelling: ovisc (sohd 
line), Ivisc (dashed line), csf (dashed-dotted), csfc (dotted line). Upper and 
lower panels refer to LM and HM samples, respectively. 



re-simulations, infact, have been chosen to be relaxed and isolated, 
while the LM clusters have more varied dynamical configurations. 
As also pointed out bv lRasia et aLl ( l2004l) . the density profiles in the 
external regions are on average steeper for isolated systems than for 
objects having more general environmental situations. This prop- 
erty, then, reflects in our SZ profiles. Specifically, profiles of both 
subsamples have a logarithmic slope going to as R/Rvir 0. 
The slope of the LM sample steepens at R/R^ir — 0.04 — 0.07 
and then remains nearly constant up to Rvh- Conversely, the slope 
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log(hM/Mo ) 



15.5 



Figure 4. Mass dependence of the fitting parameters rp (upper panel), K 
(central panel) and b (lower panel). Different symbols refer to simulations 
with different physical processes included, as indicated in the upper panel. 



of the profiles in the HM sample presents a first steepening at 
R/R^ir = 0.05 - 0.09 and a second one at R/R^ir ~ 0.6. To 
reproduce this feature we had to add a negative constant {K) to the 
fitting relations. 

This result suggests that the fitting parameters Vp, b and K 
might depend on the cluster mass. To investigate whether this is 
true, we fitted relations like Eq.(l5]( to all individual SZ profiles, and 
plotted in Fig. |4] the resulting parameters vs. cluster mass. Since 
the mean profiles depend only mildly on the cluster physics, we 
prefer to show the data coming from the different re-simulations 
together, thus increasing our statistics. Plotting the parameters Vp, 
K and h against the (logarithm of the) mass we do not find any 
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10-3 



10-5 



10- 



lO-'^ 



10-5 



10-6 



Ivisc vs ovisc 




csf vs ovisc 




csfc vs ovisc 



10-6 10-5 10-4 10-3 

yo 

Figure 5. Comparison of the values yo obtained for a given cluster in the 
ovisc simulation (horizontal axis) with the coiresponding values obtained 
in the other simulations (vertical axes): Ivisc (top panel), csf (central panel) 
and csfc (bottom panel). Different symbols refer to clusters observed at dif- 
ferent redshifts: 2 = (crosses), z = 0.2 (stars), z = 0.4 (squares), 
z = 0.7 (diamonds) and z = 1 (triangles). The solid line corresponds to 
the bisecting Une. 



statistically significant trend. However, we notice that, as suggested 
by the comparison of the mean profiles, higher masses correspond 
on average to higher values for the scale radius rp and to steeper 
profiles at i?/i?vir 1, as shown by the higher (absolute) value 
assumed by K. 



5 SZ SCALING RELATIONS 

If gravity is the only physical process in action during their for- 
mation, galaxy clusters can be modelled using th e so-called self- 
similar model i lKaisej|l986l :l Navarro et al. which predicts 
simple scaling relations between the main cluster properties (mass, 
temperature, etc.). In particular, in terms of SZ and X-ray observ- 
ables, the self-similar model suggests the following relations: 



yo 

yo oc T'^''E{z) , 

AS oc r''/^£-'(2) 



(6) 
(7) 
(8) 



In the previous equations yo is the central value for the Comp- 
tonization parameter y, AS is the SZ flux integrated out to a given 
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Figure 6. As Fig|5] but for AS. 



radius, T is the ICM temperature, Lx is the X-ray luminosity. The 
factor E{z), representing the ratio between the Hubble constant at 
redshift z and its present value, is given, for flat cosmologies, by 



csfc vs ovisc 




E^{z) = nomii + zy + nA 



(9) 



and allows rescaling all the data to the same redshift (usually z — 
0). 

The SZ self-similar scaling relations ha ve been largely in- 



vestigated with semi-analytic models (see, e.g., Cavaliere & Mencj 



m numer- 



I2OO1I : I Verde et alj |2002| : iMcCarthv et al] |2003|) and 

ical hydrodynamic simulations (see , e.g., Ida Silva et alj |2004| : 
lOiaferio et al.ll2005l : iMotl et aljlioo^ . since deviations from them 
would indicate that non-gravitational processes, like gas cooling 
and energy feedback, play an important role in gas-dynamics. It is 
worth noticing that strong support in favour of a non-negligible role 
of the non-gravitational processes already comes fr om the study of 
self-s i milar relations in the X-ray band (see, e . g. ICavaliere et all 
'1998"; 'Allen & Fabian" 1998"; 'Mai-kevitch' 1998; Xue & Wu'l200(:J; 
Borgani et al. 2004; Arnaud et al. 2005; Vikhlinin et al. 2006, and 
references therein). 

We studied the relations in Eqs.l|6j-l[8j by analysing the same 
cluster sample of the previous section, but considering each cluster 
at five different redshifts: z = 1, 2 = 0.7, z = 0.4, z = 0.2, z = 0. 
This corresponds to a time separations ~ 2 Gyr, long enough to 
allow changes in the cluster dynamical status, and thus to make 
successive data points from the same cluster roughly independent 
on each other. 

To construct our empirical relations, we used the intrinsic pa- 




Figure 7. Scaling relation yo — Lx for all physics. Solid and dotted lines refer to our best-fit relation and to the best fit obtained by imposing the self-similar 
slope, respectively. Different symbols refer to clusters at different redshifts, as in Fig. [5] 



Table 4. Best-fitting parameters for the different SZ scaling relations (j/o ~ 
Lx, yo — T and AS — T) obtained from our hydrodynamic simulations 
with different physical modelling of ICM. X-ray luminosity in the [0.1-10 
keV] band Lx, mass- weighted temperature T and integrated SZ flux AS 
are given in units of lO*** erg/s, keV and mJy/Mpc^, respectively. 



Physics 




yo - Lx [10 


*■* erg/s] 


ovisc 
Ivisc 
csf 
csfc 


log{j/o) = 
log{j/o) = 
log(?/o) = 
log(?/o) = 


(-4.77 ±0.01) -f 
(-4.76 ± 0.01) + 
(-4.63 ± 0.01) + 
(-4.59 ± 0.01) + 


(0.78 ± 0.02) log(Ljf) 
(0.77 ± 0.02) log(Lx) 
(0.79 ± 0.02) log(Lx) 
(0.81 ± 0.03) log(Lx) 


Physics 




yo - T [keV] 


ovisc 
Ivisc 
csf 
csfc 


log(j/o) = 

log(yo) = 

log{j/o) = 
log{j/o) = 


= (-4.70 ±0.03)- 
= (-4.80 ±0.03)- 
= (-5.11 ±0.03) - 
= (-5.16 ±0.01) - 


1- (1.54 ±0.07) log(T) 
1- (1.42 ± 0.06) log(T) 
1- (1.71 ±0.05) log(T) 
1- (1.75 ±0.03) log(T) 


Physics 




AS [mJy/Mpc2] 


-T [keV] 


ovisc 
Ivisc 
csf 
csfc 


log(AS) 
log(AS) 
log(AS) 
log(AS) 


= (1.68 ± 0.01) -( 
= (1.71 ±0.02) -( 
= (1.51 ± 0.01) -( 
= (1.53 ± 0.01) -( 


- (2.40 ± 0.01) log(T) 

- (2.37 ± 0.02) log(T) 

- (2.50 ± 0.01) log(T) 

- (2.46 ± 0.02) log(T) 



rameters extracted from our SZ maps. The X-ray luminosity in the 
[0.1-10 keV] band Lx and the integrated SZ flux AS were com- 
puted within i?5oo. To estimate the temperature we preferred to 
adopt a mass-weighted estimator: this quantity is in fact related to 
the energetics involved in the process of structure formation. As 
shown in earlier papers | Mathiesen & Evrard 200l]: lMazzotta et^ 
l2004lGardirir et al.ll2004inRasia et alil2005i) . the alternative of us- 
ing an emission-weighted temperature - originally proposed to 
make results from hydrodynamical simulations directly compara- 
ble to observational spectroscopic measurements - actually intro- 
duces systematic biases when the ICM structure is thermally com- 
plex. 

First we wish to investigate how the SZ signal is modified by 
the inclusion of different physical processes. In Fig. [5] we make a 
one-to-one comparison of the value j/o obtained for each cluster 
at a given redshift in the ovisc runs versus those obtained in runs 
including more physics. A similar plot, but for AS, is presented in 
Fig|6] Different symbols refer to different redshifts, as indicated in 
the figure captions. We find that the values measured in the csfc and 
csf runs are systematically lower than for ovisc, by a factor ~ 60 
per cent for yo and of 25 per cent for AS*. On the other hand, we 
tested that the plotted data can be fitted with a logarithmic slope 
which is always compatible with unity, showing that the slope of 
the relations is robust against changes in the modelled physics. 

Let us now turn to the SZ scaling relations. The slope and 
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Figure 8. As Fig|7] but for the scaling relation ya — T. 



the noimalisation of our best fit relations, with the corresponding 
errors, have been obtained by applying to each set of parameters 
(x, y) a "robust" least absolute deviation method. Note that here 
X and y actually express a log-log relation. First, we fitted the data 
assuming a linear model 1/ = Ai+Bix; we thenrepeated the linear 
fitting swapping x and y, and obtained the new parameters {A2 and 
B2)- Finally, the slope and intercept [A, B) and the corresponding 
errors (AA, A_B) were calculated as follows: 

A=(Ai+A2)/2 , AA^ jyli - A2i/2 (10) 
B = {Bi+B2)/2 , AB = |Bi - B2 1/2. (11) 

Our best fit relations, reported in Table |4l are also plotted as 
solid lines in Figs. [Till] and |9] for the scaling relations between yo 
and Lx, yo and T, and AS and T, respectively. For comparison, 
in each plot the dotted line is the best fit when the self-similar value 
for the slope is imposed. We first notice that data belonging to dif- 
ferent redshifts and corrected by the appropriate factor E(z) define 
a tight relation: this indicates that the evolution in redshift is in 
good agreement with predictions of the self-similar model, at least 
up to 2; = 1. Since our scaling relations were obtained using a suf- 
ficiently large number of points, the errors associated to the best 
fitting slopes and normalizations are relatively small: this increases 
the statistical significance of our results. 

The empirical yo — Lx relation for all physics is consistently 
slightly steeper than the self-similar slope, 0.75; the maximum dif- 



ference is ~ 0.05, coming from simulations which include radiative 
processes. The recovered yo — T relation is shallower than the self- 
similar relation for Ivisc, while for all other physics it is steeper by 
0.1 — 0.2. Finally, the AS — T relation is very close to the self- 
similar one (having slope 2.5) for the csf runs, but shallower (by 
0. 05-0.1) for the other ICM models; this conflicts with the work 
bv lNagail ([2006), who finds indications (though with no statistical 
significance) of a slope steeper than the self-similar one. 



In general, comparing the values in Table |4] we find that the 
slope of the relations depen ds very little on the cluster physics. 
Similar results were found by|N agail j2006h . who reports that intro- 
ducing radiative cooling lowers the normalization of the AS — T 
relation, but has no significant impact on the slope. In particular, he 
finds a normalization lower by 34 per cent, higher than our result. 
Moreover, we notice that also the Ivisc clusters exhibit a normaliza- 
tion lower than the ovisc ones: part of the pressure is in fact missing, 
because it was converted into turbulent motions. Note that this pro- 
cess is completely different from the pressure reduction in the runs 
csf and csfc, where cooling removes gas from the hot phase in the 
cluster centre. Of course, it is possible that in nature both processes 
are acting at the same time. Interestingly, we find that at the cen- 
tre the signal reduction produced by turbulence and by cooling are 
comparable in high-mass clusters. 
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Figure 9. As Fig|2] but for the scaling relation A5 — T. 



6 MOCK OBSERVATIONS 

It is now important to verify if the small signatures produced by the 
different physical processes included in the hydrodynamic simu- 
lations discussed above can be detected with a significant signal- 
to-noise ratio using current SZ instruments. In order to do that, 
we produced mock observations of our clusters. Amo ng the many 
dedicated instruments which are planned (see, e .g., iMohr et al.l 
l2002l :lLol [20o3 : lKosowskvll2003l : lRuhl et al.ll2004h . we decided t o 
consider the Arcminute MicroKelvin Imager (AMI; Jones 2002|), 
which is in its commissioning phase in Cambridge and has 

already obtained a first detection of a SZ d ecrement pointing to- 
wards the cluster A1914 tearker et id]|200 ^. More in detail, AMI 
is a compact array of 10 small (3.7m) dishes packed together and 
having a low noise, high-bandwidth back-end system. A comple- 
mentary array of 8 much larger (13m) dishes will allow the detec- 
tion and subtraction of unresolved radio sources. Both arrays op- 
erate at frequencies of 12-18 GHz. The field of view is 21 x 21 
arcmin. 

Obtaining a mock observation with the AMI characteristics re- 
quires simulating the response of the interferometric array to some 
signal on the sky. To do this, we Fourier transformed the image of 
the simulated cluster, and sampled in the u, v plane according to 



see, for more detail, http://www.mrao.cam.ac.uk/telescopes/ami/ 
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Figure 10. Typical AMI n, v sampling (courtesy of K. Grainge) adopted 
for our mock observations. It coiresponds to a field at dec 50° observed at 
HA=[-4,+4] hours. Axes labels are in units of the wavelength. 



the array configuration. The it, v sampling adopted for our obser- 
vations, and shown in Figure[To] was kindly provided by the AMI 
collaboration: it corresponds to a field at dec 50° and HA [-4, +4] 
hours. 



SZ profiles and scaling relations: modelling effects and observational biases 1 1 



We then converted our SZ maps from the Comptonization pa- 
rameter J/ to a measured flux decrement (given in mjy). We subse- 
quently performed the samphng in the u, v plane and added noise 
at the appropriate level. Following Kneissl et al. (2001), the noise 
resulting in the image is AS'rms = 20 mJy s~^^^ over the field of 
view here considered. In order to reduce noise and spurious struc- 
tures due to the sampling, we smoothed the map and applied the 
CLEAN deconvolution algorithm in order to detect the cluster sig- 
nal. The resulting map, calibrated in flux for the instrumental effi- 
ciency, is what we call "observed" cluster. 

Notice that in this analysis we do not consider any bias due 
to CMB fluctuations. The spectrum of primary CMB anisotropies 
may possess significant power on the largest scales probed by AMI; 
however, on scales of few arcmin the SZ signal clearly dominates. 
As we will see, in our case we have good resolution (~ 2 arcmin), 
and we consider only bright objects, so including the effect of the 
CMB should not really change our results. 

For this observational-like approach we considered the sim- 
ulated clusters of the HM sample only: in fact, in order to get a 
signal-to-noise ratio (S/N) good enough to allow robust extrac- 
tion of radial profiles, the exposure time for small-mass objects 
would have to be too high. Moreover, we considered the objects 
at 2 — 0.2: at this redshift our clusters cover a significant part of 
the field, allowing the instrument to resolve even small structures. 
The integration time of each simulated observations is t — 34/i, 
translating into detections at 10 to 20cr significance. As an exam- 
ple, in Fig[TT]we compare one map of the SZ decrement extracted 
directly from the hydrodynamic simulation to the corresponding 
one obtained at the end of the AMI observational procedure. The 
figure refers to the ovisc run of cluster g72a, which undergoes an 
important merging event at z ^ 0.2. It is evident that in the ob- 
served map many details are completely lost and isoflux contours 
appear more regular. Moreover, the value of SZ decrement at the 
centre is reduced, as a result of the reduced resolution. 

After preparing the individual observations, we extracted the 
mean profiles as described in the previous sections. The results for 
clusters simulated including different ICM modelling are shown 
by the thick lines in Fig. [T2] In the same plot, for comparison, we 
also display as thin lines the mean profiles directly extracted from 
the hydrodynamic simulations. The ripples in the observed profiles 
are due to the fact that our simplified removal procedure could not 
completely cancel spurious structure introduced by the sampling. 

As expected, the observational procedure makes it impossible 
to appreciate differences between SZ profiles produced consider- 
ing different physical processes. The observed profile is modified 
by the instrumental beam, so the physics mainly affects the total 
intensity of the emission. However, this quantity also depends on 
the cluster mass, and breaking the degeneracy is difficult. 

In Table[5]we discuss the robustness of the parameters j/o and 
AS used in the scaling relations, against the observational process. 
The central Comptonization parameter yo is clearly suppressed by 
the degraded angular resolution; the amount of suppression de- 
pends both on the physics included and on the shape of the specific 
cluster, ranging from a factor of 1.5 to almost a factor of 2.6 for 
the ovisc run of cluster g72a. Conversely, the parameter AS is less 
prone to observational effects; moreover, the amount of this effect 
is very similar for all clusters and for all physical processes, be- 
ing always ~20 per cent. We can thus say that, with respect to the 
observational process, the parameter AS is much more stable than 
j/o- As a consequence, the results we found for the AS — T relation 
should hold even if a properly observed sample of clusters is used. 



Table 5. Ratio between the SZ quantities (yo and AS) extracted directly 
from simulated clusters and from their mock observations (i/oobs ™d 
A5obs)- Columns refer to the different galaxy clusters belonging to HM 
sample, while rows show the ratios for simulations including different sets 
of physical processes. 
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Figure 12. The mean logarithmic SZ profiles obtained from mock ob- 
servations of simulated galaxy clusters are shown as thick lines. Different 
hues refer to simulations including different sets of physical processes: line 
types ai'e as in Fig[3] For compaiison, the thin lines refer to the coirespond- 
ing mean logarithmic SZ profiles directly extracted from the hydrodynamic 
simulations. Only the HM sample aX z = 0.2 is considered here. 



7 CONCLUSIONS 

We used a sample of simulated galaxy clusters with masses rang- 
ing from 10^"* to about 2 x 10^^/i~^Mq, re-simulated using dif- 
ferent sets of physical processes (radiative cooling, star formation, 
feedback and thermal conduction) to investigate the dependence of 
the SZ emission on the ICM physics. Starting from the 3D pres- 
sure information, we obtained for each physics fitting relations 
to the SZ radial profiles. We compared the scaling relations link- 
ing cluster central and total SZ flux, and cluster X-ray luminosity 
and (mass-weighted) temperature, with those predicted by the self- 
similar model. We then performed mock observations of the most 
massive systems in our sample, in order to check whether our find- 
ings survive the observational procedure. We summarize our main 
results below. 

• The shape of the SZ cluster profile strongly depends on the 
cluster mass, as a result of the fact that our high mass clusters, at 
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Figure 11. Maps for the SZ decrement for the ovisc simulation of cluster g72a. The original map extracted from the hydrodynamic simulation, and the same 
map in the simulated observation (t = 34/i) which assumes the AMI interferometric response, are shown in the left and right panel, respectively. The side of 
each map corresponds to 16 arcmin. The colour scale is shown at the bottom of each panel. 



variance with the low mass ones, are all isolated and virialized. 
In particular, while the 3D (logarithmic) pressure profiles of low- 
mass clusters (LM sample) can be described by a function which 
is nearly flat for i? ^ and steepening for R/Rvir ~ 0.05, the 
high-mass clusters (HM sample) exhibit an additional steepening 
of the profiles at R/R^ir ~ 0.6. 

• The details of the ICM physics have little impact on the SZ 
profiles; the largest differences are found near the central regions 
of low-mass systems. 

• The scaling relations between cluster SZ emission and cluster 
X-ray luminosity and temperature evolve with redshift according 
to the self-simil ar model out to z — 1; this result is consistent wit h 
previous works jda Silva et al .l2004lMotl et al.l2005LlNa"gaj2006h . 

• The slopes of the scaling relations for all physics are consis- 
tent with the predictions of the self-similar model, even if small 
discrepancies can still be found: in particular, our empirical rela- 
tions 2/0 — Lx and yo ~ T are generally slightly steeper than the 
self-similar ones, while the AS — T relation is generally slightly 
shallower. 

• The ICM physics has little effect on the slopes of the scal- 
ing relations , in agreement with previous studies jWhite et alj2002l : 
lNagaill2006h . 

• The main effect introduced in the scaling relations by radiative 
cooling is a lower normalization. If we compare our simulations in- 
cluding cooling with the non-radiative runs with standard viscosity 
(ovisc), we find a suppression of 50 per cent in the central intensity 
2/0 and of 25 per cent in the integrated flux AS. If we compare the 
cooling plus star-formation runs with the non-radiative ones adopt- 
ing modified artificial viscosity (Ivisc), these reductions are of the 
order of 20 — 30 per cent. A similar norma l ization lowering for the 
AS ~ T relation was also found bv iNagail ( l2006h . 

• Our mock observations suggest that arcminute resolution data 
from the current generation of SZE imaging experiments will not 
be a sensitive probe of the different physics that is ongoing in the 



intracluster medium. At arcminute scale angular resolution these 
effects tend to only have observable consequences on the total clus- 
ter flux, which leads to a degeneracy between physics and cluster 
mass that will be difficult to break 

• The integrated flux AS is found to be stable against the obser- 
vational process, and thus should be preferred to j/o, which instead 
exhibits large variations after performing the mock observations. 
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